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Abstract -We describe a scheme with analytic result that allows to generate steady-state entan- 
glement for two atoms over a dissipative bosonic medium. The resonant coupling between the 
mediating bosonic mode and cavity modes produces three collective atomic decay channels. This 
dissipative dynamics, together with the unitary process induced by classical microwave fields, drives 
the two atoms to the symmetric or asymmetric entangled steady state conditional upon the choice 
of the phases of the microwave fields. The effects on the steady-state entanglement of off-resonance 
£\J ' mediating bosonic modes are analyzed. The entanglement can be obtained with high fidelity re- 

gardless of the initial state and there is a linear relation in the scaling of the fidelity with the 
\Q • cooperativity parameter. The fidelity is insensitive to the fluctuation of the Rabi frequencies of the 

IQ I classical driving fields. 

^ ' 

in 
o 

(N , 

Introduction. — Quantum entanglement plays a crucial role in performing quantum 
information processing [IHl], and it is the basic requirement for quantum communication 
to transfer quantum state between distant locations. To construct quantum networks, one 
; I \ needs to generate spatially separate qubits, store them for sufficiently long time, and perform 

logic operations between them. Therefore, the generation of entanglement between different 
nodes becomes an important task [5]. The optical absorption and other channel noise 
inevitably bring the detrimental effects to the entanglement between distant nodes, which 
normally decreases exponentially with the length of the connecting medium. The quantum 
communication difficulty cannot be solved just based on unitary dynamics [Sttldj . 

To overcome the problem associated with the exponential fidelity decay of distributed 
entanglement, the concept of dissipative dynamics can be used, which utilizes the dissipa- 
tion as a resource for quantum state engineering [?. ?. ?. ?. [T51[TSH^51[?7H551l37j . Recently, 
experimental realization of dissipative state preparation has also been reported |38j. Here, 
we describe a scheme with analytic result that allows to generate steady-state entanglement 
for two atoms over a dissipative bosonic medium. Compared with the previous scheme we 
proposed in Ref. |18j , the steady entanglement can be produced over a longer distance via a 
dissipative mediating mode, and is very useful for testing quantum nonlocality [39], quantum 
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Fig. 1: (Color online) Experimental setup and level diagram for dissipative preparation of long- 
distance entanglement between two atoms coupled to two cavities respectively, which are connected 
by a multi-mode dissipative medium. 

secret sharing |40j and quantum communication |41j . 

Due to the coherent cavity-medium coupling, the system becomes mathematically equiv- 
alent to the case that two distant atoms collectively interact with three nondegencratc delo- 
calizcd field modes, one mode inducing the asymmetric collective atomic decay channel and 
the other two inducing the symmetric collective atomic decay channels. The decay chan- 
nels induced by three delocalized field modes are analytically resolved with the dissipative 
dynamics. Analytical result indicates that the distributed entanglement can be obtained 
with high fidelity requiring neither the unitary feedback control nor the photon detection. 
The steady-state entanglement arises from the competition between the collective atomic 
decays and the unitary evolution induced by the classical microwave fields, which act as a 
push button to start the dynamical process. We find that the scaling of fidelity F with the 
cooperativity parameter C is linear. Besides, the effect on the steady-state entanglement 
of dispersive mediating modes analyzed. We find that the couplings between the dispersive 
mediating modes and cavity modes lead to the frequency shifts of the normal delocalized 
bosonic modes, when one of these modes is resonant with the classical laser field, the Raman 
transitions between two atomic ground states dominate the atom-field coupling dynamics 
and destroy the process of the steady-state entanglement preparation. The steady-state en- 
tanglement is faithfully prepared when the frequency space of the mediating bosonic modes 
is large enough, showing that our scheme can work well even when many mediating bosonic 
modes are involved. 

The theoretical model. — As shown in Fig. 1, two three-level atoms are individually 
trapped in two single-mode cavities linked by a dissipative bosonic medium (a third cavity 
or a waveguide). Each coherent driving D, pumps the ground state |0) to the excited state 
1 2) with the detuning A. Atom-cavity interaction g drives the transition |1) -B- |2) with the 
detuning A — S, and two ground states are coupled by a resonant microwave Qm in each 
cavity. 

We assume two cavity modes are coupled to the mediating mode. Then the Hamiltonian 
of the whole system in the interaction picture can be written as Hi = Hq + H g + V+ + V- : 
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H g = ^(|l)i<0|+e- w *|l) a <0|) + iT.c, (2) 
V + = ^(I2>x<0| -I- 12> 2 <0|), (3) 
V- = ^(|0) 1 (2| + |0) 2 (2|), (4) 

where a 2 : (i = 1,2) is the annihilation operator for the ith cavity field mode, b n is the 
annihilation operator for the nth mediating mode, A„ denotes the frequency difference 
between the nth mediating mode and the cavity mode. v n denotes the coupling strength 
between the ith cavity mode and the nth mediating mode, g is the atom-cavity coupling 
constant, Q and Om represent the classical laser driving strength and the microwave driving 
strength, respectively. 9m is the phase difference between the two microwave fields applied 
to atoms 1 and 2. The cavity- mediating couplings are assumed to be the same for all the 
mediating modes, i.e., v n = v. 

To analyze the main dissipation mechanism easily, we first focus on the case that only 
one mediating mode (i.e., the first mediating mode b\) resonantly interacts with the cavity 
mode and the other mediating modes with large detuning can be neglected in Hq. Let 
us define the symmetric steady-state |T) = (|10) + |01}) / a/2, the asymmetric steady-state 
IS*) = (| 10) — |01))/-\/2, and introduce three normal delocalized bosonic modes as C\ = 
(ai — a 2 )/v / 2, c 2 = (ai + a 2 + \/2&i)/2, and C3 = (ai + a 2 — \f2b\)/2. In terms of the new 
operators, the Hamiltonian Ho becomes 

H Q = A(|2) 1 (2| + |2) 2 (2|) 

9|2)i(l|(ic 2 + ic 3 + ^ Cl ) 

+ g\2) 2 (l\( 1 -c 2 + l -c i -^ Cl ) + H.c. 

+Sc\ Cl + {5 + \/2V)4c2 + {S - v / 2^)4c 3 . (5) 

Here, Hq describes the asymmetric coupling for two atoms to the delocalized field mode 
Ci, and the symmetric couplings to c 2 and C3. Because of the cavity- medium coupling, 
all the delocalized field modes arc nondegenerate and each causes a collective dissipation 
channel, (i = 1,2) and K3 represent the photon decay rate of cavity i and the medium, 
respectively. 70 and 71 denote the atomic spontaneous emission rates for |2) — > |0) and 
|2) — > |1), respectively. For simplicity we here set «i = k 2 = k 3 = k and 70 = 71 = 7/2. 
Then the corresponding Lindblad operators associated with the photon decay and atomic 
spontaneous emission can be expressed as L Kl = ^fti c%, L K2 = ^/k c 2 , L K3 = ^/k C3, U 11 = 
V7? |0)i(2|, = |0) 2 (2|, = ^ |l)i(2|, = VtT |l)a<2|. 

Under the condition that the classical laser fields are weak, we can apply the method 
based on the second-order perturbation theory to tailor the effective dissipation process 
[15H17] . Then the dynamics of our distributed cavity QED system is governed by the 
effective Hamiltonian H e ff and effective Lindblad operator L^^: 

Heff = -~V-[Hj? H + (Hj H )i]V + + H a , (6) 

L eff = LXH NH V +> ( 7 ) 

where H^ H is the inverse of the non-Hcrmitian Hamiltonian Hnh = Hq — ^^2 X (L X Y L x , 
and its elements determine the strength of effective evolution process. The resulting effective 
master equation in Lindblad form becomes 

P = i[p,Heff]+J2i L effP( L effV 

X 
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-\[{Ll f ^L% n p + p{Ll f ^Ll f} }}, (8) 
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Hgfj = -ReltfRi] |00) (00| - Re[—R 2 ]\S) (S\ 

-Re[^R 3 ]\T)(T\ + H g , (9) 
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where Re[ ] denotes the real part of the argument, 



g e = gCl,8 =8 -in/2, A = A — ry/2, 
5'{8' 2 -2v 2 ) 



Ri 
R 2 

Ra 



A'8'(S' 2 -2v 2 )-g 2 (S' 2 -v 2 y 

A'S'(S' 2 -2v 2 )- g 2 S' 2 
(g 2 - S'A')(S' 2 A' - S'g 2 + 2AV) ' 

AV(,5 ,2 -2i/ 2 )-. 9 2 ((5 ,2 -2z/ 2 ) 
(g 2 - S'A')(S' 2 A' - S'g 2 + 2AV) ' 
A = A8{8 2 -2v 2 )-g 2 {5 2 -v 2 ), 
B = {5/2-v 2 ){Ak + 1 5) + 8k{A5 - g 2 ), 
d = g 2 -A8, 
Dt = (Are + <5 7 )/2, 
C 2 = g 2 S- A{5 2 -2v 2 ), 

D 2 = n{A8-g 2 /2)+ 1 {5 2 -2v 2 )/2. (13) 

Since f2 is very small, we can neglect those terms containing the effective shifts 0(f2 2 ) in Eq. 
(9), so that H e ff ~ H g . There are three primarily effective decay channels characterized 
by L" 1 , ig 2 , ig 3 through three dclocalized bosonic modes c\, c 2 , C3, respectively. It is the 
mediating mode that links the two separate field modes and lifts the degeneracy of three 
dclocalized field modes, leading to three independent decay channels. We assume re, 7 <C g, 
5, v, A, so that the minor terms with higher order than re 2 , j 2 and K7 in the denominators 
of Lg 1 , Lg 2 and L* 3 have been omitted. 

As shown in Figs. 2(a) and 2(b), L" 1 indicates the effective decay from 1 00) to \S) at a 
rate re* and from IS) to |11) at a rate re 2 ^, which is caused by c\ mode that only contains 
the contribution of cavity mode. Lg 2 (L* 3 ) indicates the effective decay from |00) to |T) at a 
rate re* 2 (re* 3 ) and from |T) to |11) at a rate k 2 x (re 2 J, which is caused by c 2 (03) mode that 
contains the contributions of both the mediating mode and cavity modes. Three effective 
collective decays happen simultaneously. The decay rates re 1 ^ (re 2 J, re* 2 (re 2 2 ) and re* 3 (re 2 3 ) 
equal to the square of the first (second) coefficient in the right hand side of Eq. (10), Eq. 
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Fig. 2: (Color online) Effective two qubit system for coherent and dissipative processes in the ground 
state basis | |00), \S), \T), 11) }. Effective dissipative rates induced by the field decay from 1 00) 
to \S) and |T) are and (nl 2 + kJ 3 ), and from |5) and \T) to |11) are an d (k 2 2 + k% 3 ). 7g 2 , 
7g 4 , 7y 2 , and 7^,' 4 are the effective dissipative rates due to atomic spontaneous emission, (a) 8 = 
tt; (b) 6 = 0. 



(11), and Eq. (12), respectively. To prepare the desired steady-state entanglement \S) (|T)), 
it is necessary to suppress the effective decays from \S) to |11) (from \T) to |11)), i.e., the 
condition n\ x >• n 2 Ci (k^ 2 >• k 2 2 and k\ 3 ^> k^) should be satisfied. Two microwave fields 
cause the coherent shuffling of the three states 1 00} , \T) (\S)) and |11) for 6 m = (tt). As 
a result, the state \S) (\T)) is the unique stationary point of the system. 

The fidelity. — Using Eq. (7), we obtain the effective Lindblad operators that drive 
the population out of the target state \S) 

LZ]s = ^ s = yfff\T)(S\, (14) 
LT, S = I% S = ^\U)(S\, (15) 
and those driving the population out of the target state |T) 

L^T = L^ T = yf^\S)(T\, (16) 

L?,t = L^ T = ^\ll){T\, (17) 

where 7e = 7 ft 2 ^ 4 / [ A(8 2 -u 2 ) + 5g 2 ] 2 , 7^ = 7 ^ 2 = 7e /32 and 7 §' 4 = 7 ^' 4 = 7e /16. The 
evolutions of populations of different state components governed by the full master equation 
are plotted in Fig. 3, which shows that the steady-state entanglement \S) or \T) can be 
obtained with high fidelity depending upon the choice of the phase difference 9m between 
the two microwave fields. 

To evaluate the fidelity of the steady-state \S) or |T), we apply the rate equation as 
follows 

Ps = <^oo-«+7s 2 +7s 4 )^, (18) 
Pt = «+<)Poo 

-«+<+7^ 2 + 7^ 4 )Pt, (19) 
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Fig. 3: (Color online) The populations of states 15), |T), 1 00) , and |11) versus the dimensionless 
parameter gt for a random initial state by solving the full master equation. The curves are both 
plotted for an appropriate set of parameters C = 150, 7 = 2k, k — 0.0577ij, 7 = 0.1154<?, O = 
O.O63, fi M = 0.0138s, A = 1.30, v = 0.4528ff and S = 0.2875p. (a) 6m = 0; (b) 6 M = ir. 



where P v is the probability to be in the state \<p). When the whole system reaches the 
steady entanglement \ip), we assume P v ~ 1 and the probability in each of the other three 
states is nearly P 00 , then 1 - F\ B) « 3P o (B = S, T). Here, F\ B) = \ (B\p\B) \ is the fidelity 
of state B, and 

f 3 ge 2 ^ 9 7e . (6 2 - 2^fg 2 e n 
[ 8(C 2 + Dl) 32 1/1 4B 2 J ' 

(20) 

3(£ 2 + 2^)g e 2 K 9 7e * 2 (* 2 + 2i/ 2 W 
[ 8(C 2 + L> 2 ) 32 1/1 AB 2 J ' 



1-F, S) « 
1-F m » 
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Fig. 4: (Color online) The fidelity for steady-state entanglements \S) (red solid line) and \T) (black 
solid line) vs the cooperativity parameter C, and the mathematical fitting curves for maximizing 
F\s) (blue dot line) and F\t) (green dot dash line). 



(21) 

We plot the fidelity for entangled steady-states with the cooperativity parameter C and gain 
a simplified formula through mathematical fitting based on the least squares algorithm in 
Fig. 4, and find out the actual constants for maximizing the fidelity as follows 

l-ifa ~ 14.5C7" 1 , (22) 
l-F\ T) ~ 27.5C" 1 . (23) 

The fidelities of the steady-state versus the fluctuations of coherent driving and mi- 
crowave fijvf are plotted with full master equation in Figs. 5(a) and 5(b), and the results 
shows that the fidelity F\$) and F\t) remain higher than 91% even when the relative errors 
in fl and SIm reach 20%. 




6QJQ. 6Q.IQ. 

(a) (b) 

Fig. 5: (Color online) Fidelity of the steady states versus the fluctuations of Q and Qm. d£l and 
dQ.M are the corresponding deviations of Q and Qm, respectively. Both plots are calculated with 
full master equation at time t = 9000. (a) F\g\ vs dil/Q and dflM /Qm; (b) Fit) vs d£l/Q and 
dQ.M /SIm ■ 
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Influences of the dispersive mediating bosonic modes. In the above deriva- 
tions, we assume that only one resonant mediating mode interacts with the cavity modes. 
We now consider the effects on the steady-state entanglement of other existing dispersive 
bosonic modes. Fig. 6 shows influences of the number ./V of the mediating modes and the 
frequency spacing A x on the fidelity \F)s- 

F\t) is influenced severely by the dispersive mediating modes and it is difficult to obtain 
the steady-state \T) when N becomes large. Therefore, we just discuss the state IS*) as 
follows. We first consider the case that two mediating modes are involved in Fig. 6(a), 
where we assume the first mediating mode is resonant with the cavity modes. There is a dip 
when A x = — 0.64g. This can be explained as follows. The coupling between the mediating 
modes and cavity modes leads to frequency shifts of the delocalized field modes. When 
one of the delocalized modes, together with the classical laser field happens to induce the 
resonant Raman transition between two ground states (|0) and |1)), as shown in Fig. 7, the 
Raman process, together with cavity decay, leads to effective decay |0) — > |1) and destroys 
the entanglement. 

The fidelity versus frequency spacing with three (five) mediating modes is plotted in Fig. 
6(b) ( Fig. 6(c) ), showing that two valleys appear when A x = — 0.54(7 (— 0.58g) and = 
0.54g (0.58<?). The appearance of two symmetric valleys in Fig. 6(b) ( Fig. 6(c) ) is due to 
the fact that two (four) equally spaced mediating modes are symmetrically distributed in 
two sides of the resonant mode. Especially, two small valleys appear when |A X | = 0.20g in 
Fig. 6(c), where two of those delocalized field modes are degenerate. 

Numerical simulations show that the fidelity of the steady-state entanglement is free 
from the effects caused by dispersive mediating modes when A x > 6. Strong coupling for 
photons between distant nanocavities mediated by a waveguide at room temperature has 
been reported and the system dynamic control can be realized |42j . The cooperativity 
parameter C about 100 has been realized [43]. Based on these experimental techniques, the 
present scheme is feasible. 

Conclusions. — In summary, we have proposed a dissipation based scheme by which 
two atoms trapped in two distant cavities can be driven to steady-state entanglement. The 
competition based on the unitary dynamics and the dissipative dynamics leads to the en- 
tangled steady-state. Our results show that the dissipative mediating bosonic mode can be 
used as an entanglement catalyst for two distant atoms. The effects of dispersive mediating 
modes are analyzed. It is shown that the entanglement is robust against the fluctuations 
of the Rabi frequencies of the classical fields. We show that the distributed steady-state 
entanglements can be obtained with high fidelity regardless of the initial state and there is 
a linear relation in the scaling of the fidelity with the cooperativity parameter. 
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